Surface- active dust in a plasma sheath 
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Abstract 

The inhomogeneity of the plasma pressure near a conducting electrode is a cause for introducing 
the surface tension associated with the plasma-electrode interface. We evaluate the dependence of 
the surface tension on the density of the charged dust immersed in the plasma sheath. In a wide 
range of parameters, the surface tension turns out to be an increasing function of the dust density. 
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Throughout the recent decade, an admixture of charged macro-particles, electrons, ions, 
and neutral atoms, called a dusty plasma, has been the subject to some thousands of studies. 
Numerous aspects of the physics of dusty plasmas are discussed in the series of reviews 
starting with and the monograph Q]. 

Many laboratory experiments deal with relatively small number of dust grains levitating 
above a horizontal rf-powered or dc-biased electrode. Under these conditions, dust grains 
are negatively charged and the gravity force is compensated by the strong electric field in a 
plasma sheath adjacent to an electrode. Dust suspended in a plasma sheath self-organizes 
itself in various single or multi-layered structures. In particular, recent experiments [^J 
demonstrated that with growing rf power or gas pressure there evolves a void, that is, 
a dust-free region, in the center of a single dust layer. Since negatively charged grains 
levitating at the same height repulse one another via the screened Coulomb potential, there 
are no evident reasons for the void emergence. The explanation of this phenomenon given 
in [2! implies existence of an attractive interaction between grains; however, the nature of 
the interaction remains obscure. 

To avoid confusion, it should be stressed that three-dimensional voids, which have been 
observed earlier jj, 0, Q, E | > evolve in the plasma bulk. The theory of three-dimensional voids 
demands large number of grains and strong influence of the dust upon the discharge 
structure. At first sight, both assumptions fail under the conditions of the experiment |3|, 
where some hundreds of grains only have been used. 

The main purpose of this Letter is to introduce the concept of plasma surface tension. 
We demonstrate that even a small number of dust grains can modify the plasma sheath 
profile in such a way that the sheath may confine the dust in the horizontal direction. Our 
reasoning is based on the analysis of the plasma pressure in the sheath-presheath area of a 
gas discharge. We show that the horizontal electric force appearing due to the inhomogeneity 
of a dust layer is conveniently treated in terms of the plasma surface tension. It should be 
stressed that explanation of the experiment Q] is beside our purpose: there is evidently a 
number of relevant factors that are not discussed here. 

The schematic of the model implemented here is depicted in Fig.^ The plasma consisting 
of electrons and singly charged ions is described with the help the continuity equation 

dm(z)vi(z) 

-r = Vionn e {z), (1) 
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FIG. 1: Schematic of the model 
the momentum balance equation 



dm * n f vKz) +en l (z)^± = 0, (2) 



and Poisson's equation 



= 4vre(n e (z) - m(z)) - Anp d (z). (3) 

Here, Vi on is the electron impact ionization frequency and pd{z) is the charge density of 
aerosol grains. Other notation is standard. Electrons are supposed to be Boltzmannian: 
n e (z) = no exp (e<f)(z)/T e ), where no is a normalizing constant. We also assume that the 
mean free path of the ion-neutral collisions is large compared to the Debye length, which 
determines the spatial scale of the sheath. 

Since the electron and ion fluxes at the electrode (z = 0) are equal, one of the boundary 
conditions to Eqs. (PEJ) is 



/MfUoiO) -„,(()) \/^- = 0. (4) 



Far from the electrode the plasma is quasi- neutral, i.e., n e = rii = n . Generally, the 
equilibrium plasma density, no, is the eigenvalue of the problem and it is determined by the 
balance of ionization and losses of charged particles. For this reason, the bulk plasma density 
may depend on the dust charge located in the sheath because the latter modifies the electric 
field profile, which, in its turn, influences the net particle flux at the electrode. However, here 
we assume the overall size of the plasma to be large enough, so that the equilibrium density, 
n , is conditioned by the recombination in the plasma volume. Therefore, the plasma density 
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tends to a certain fixed value at z — > oo. In order to patch the quasi- neutral plasma and 
;he sheath, the equations (I1I2|) are solved assuming that n e (z) = rii(z) and pd{z) = (e.g., 
lo| ) that results in 



where t> s = ^T e /rrii is the ion sound velocity. It is assumed that at the end of the simulation 
area, z = Zq, the ion flow velocity is fixed, i.e., v^Zq) = v . The value of Vq should be small 
enough; all solutions exemplified below are obtained for vq = —0.01v s . Then, Eqs. (jl|5|) 
are used to evaluate boundary values of the density, n^Zo), the potential, 4>(z ), and its 
derivative, <f)'(z ); this results in four boundary conditions for Eqs. (PE)). The length of the 
simulation area, zo, is determined with the help of the fifth boundary condition (|4"|). 

In the framework of our model, the dust is represented as an infinitesimally thin charged 
massive layer levitating at z — Zi, that is, its charge density is 

Pd(z) = -qcrS(z - z{), (6) 

where — q (q > 0) is the charge of a single grain and a is the surface dust density. The 
equilibrium levitation height, Zi, depends on the vertical force balance. The largest forces 
acting upon dust in a plasma sheath are the electric field force and the gravity force. Since 
the z-component of the electric field, E(z) = —(f>'(z), is discontinuous at z = z\, the vertical 
force balance is written as 

^q[E(z l + 0) + E(z l -0)]=Mg, (7) 

where M is the mass of a single grain and g is the gravity acceleration (in Fig. ^ the Earth 
is at the bottom). On the other hand, the electric field discontinuity is given by 

E(zi + 0) - E(zi - 0) = -Anqa (8) 

and the potential, <fi(z), is a continuous function at z — z%. The relations ()7|8j) allows us 
to patch the numeric solutions of Eqs. (UHHJ) below (0 < z < z{) and above (z\ < z < z ) 
the layer, and to obtain the equilibrium levitation height, zi, and the net thickness of the 
sheath-presheath area, z . 
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It is well-known that the grain charge, q, in Eqs. (|7|8 1 depends on the plasma parameters, 
and its variability should be also taken into account llj . However, here we assume that 
the grain charge is fixed; the results of preliminary computer runs with the variable grain 
charge did not alter our main conclusions. We also ignore the ion wind force in Eq. (J7J) as 
long as absorption of plasma by dust grains. Simple estimations show that these processes 
are of minor importance. 

Our main objective is to analyze the plasma pressure profile in the sheath-presheath area. 
Evidently, Eqs. (001) ensure the conservation of the net linear momentum of the system. 
The tensor of the momentum flux, or the pressure tensor is a sum of the ion momentum 
flux, the electron pressure and Maxwell's stress: P Q/ g = 77^^0,1;^ + n e T e 5 a p + E 2 5 a f3/8n — 
E a Ep/Ait (a, f3 = x,y,z). In the one-dimensional model described above, the pressure tensor 
is diagonal: P zz = Pi and P xx = P yy = P tT . The longitudinal pressure is 

P,(z) = mMzMz) 2 + n e (z)T e - -^E(z) 2 . (9) 

Since within our model there are no momentum losses in the plasma, P; is constant in 
the absence of dust and it is equal to the plasma electron pressure at z — > 00, Pi = n T e . 
The transverse pressure is written as 

P tr (z)=n e (z)T e + ^-E(z) 2 . (10) 

In contrast with Eq. (0), there is no contribution of the ion momentum flux and the electric 
field pressure is of the opposite sign. 

Although the transverse pressure tends to the electron pressure in the plasma bulk, 
Ptr\z-*oo — n oT e , it is essentially nonuniform in the sheath-presheath area adjacent to the 
electrode. This circumstance allows us to define the surface tension at the plasma-electrode 
interface as 

00 

s(a) = J dz {n T e -P tr {z)). (11) 


The physical meaning of this quantity is fairly evident. Since the bulk plasma pressure is 
positive, two parts of the discharge repulse one another. As we shall demonstrate later, 
the transverse pressure ()10|) is reduced near the wall; therefore, there arises the additional 
horizontal force attracting two parts of the plasma-wall interface. The surface tension (fTTj) 



characterizes the magnitude of this force per unit length. It should be noted that usually 
surface tension is defined in terms of surface free energy. The latter is obviously meaningless 
in application to an open, nonequilibrium system like the plasma-wall boundary. However, 
due to the momentum conservation, the surface tension given by (jllj) is a well-defined 
concept. Of particular interest for our purposes is the dependence of the surface tension 
(fTTj) on the dust density. 

The utility of the notion (jllj) may be illustrated in the following way. Suppose that the 
layer density is a weakly varying function of the transverse coordinate, x, that is, da/dx ~ e, 
where e< 1. Integrating P xx component of the pressure tensor with respect to z at nearby 
points x and x + Ax we see that the horizontal force (i.e., the momentum flux) acting upon 
a small part of the layer of the width Ax is written in two identical forms: 



Comparing both parts of the identity and noting that E x ~ d/dx ~ 0(e) we conclude 
that in order to obtain the first-order term of the expansion of the horizontal force in powers 
of e it is sufficient to evaluate the surface tension ignoring the spatial inhomogeneity 
of the layer density. Therefore, f x = ds(a(x))/dx. This simple reasoning may be confirmed 
by the more tedious asymptotic expansion of the hydrodynamic equations in powers of e. 

We can also extract some useful information about the horizontal force even for the 
sharply bounded layer. Suppose that the dust layer occupies the area x < 0, as depicted in 
Fig. ^ The levitation height, zi, now depends on x, the electric field and the ion flow near 
the edge of the layer are essentially two-dimensional. However, the horizontal component of 
the force is the difference f x = s — Si, where s ,i are given by the integral (fTTj) evaluated 
across the sheath in the dust-free area, at x = xo, and in the presence of dust, x = x\ 
(Fig. Q). Assuming that both points £0,1 are far enough from the layer edge, the horizontal 
components of the ion velocity and the electric field are negligible. Therefore, the values of 
the surface tension, so,i, which determine the net force exerted upon the edge, are provided 
by the solution of the one-dimensional problem (P|SJ). 

Before discussing the results of the numeric evaluation of the surface tension, it is conve- 
nient to introduce the normalized quantities. In the following, the coordinate, z, is measured 
in the units of the electron Debye length, = ^yT e /Ane 2 no, the pressure is normalized to 
n T e , and the surface tension is normalized to X D n T e . The dimensionless surface density of 




(12) 



6 



(a) 




FIG. 2: The pressure (a) and the electric field (b) profiles in the sheath-presheath area evaluated 
for fi = 0.5, 77 = 0.1. (al) corresponds to the longitudinal pressure ©, (a2) corresponds to 
the transverse pressure ()10|) . the dashed curve (a3) corresponds to the transverse pressure in the 
absence of dust (77 = 0), (bl,b2) correspond to the electric field for rj ^ and 77 = 0. The vertical 
dashed lines in both figures show the position of the sheath-presheath interface {z s = 10.31). 
Discontinuities at the curves (al, a2, bl) are at the equilibrium position of the dust layer (zi = 3.94). 

the dust layer is r] = qa / (2enoA/j) and the normalized grain mass is \i = MgeXo/ (qT e ). Un- 
der the typical experimental conditions, fi < 1, while 77 < 0.1. All computations discussed 
below have been performed with the sufficiently small value of the ionization frequency, 
Vion/upi = 0.01, where uj pi is the ion plasma frequency. 

Examples of the pressure and the electric field profiles are depicted in Fig. |21 The net 
simulation length for this particular run is Zq = 68.5. In order to make the effect discernible 
in the figure, we use here a sufficiently large value of the normalized dust density, 77 = 0.1. 
The longitudinal pressure drop (the curve 1 in Fig. |2^i) at the charged layer is compensated 
by its weight; outside of the layer, the longitudinal pressure is constant. Contrary, the 
transverse pressure (jlOj) is strongly inhomogeneous inside the sheath (the curves 2 and 3 in 
Fig. |2k). Of importance is that the dust reduces the transverse pressure between the layer 
and the electrode. As it is readily seen from Fig.EJo, the effect is mostly due to the reduction 
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FIG. 3: Dependence of the normalized surface tension on the dust density, 77. 1 — /i = 0.8, 2 - 
H = 0.5, 3 — \x = 0.1, 4 — n = 0.02. 




FIG. 4: The derivative of the surface tension with respect to the dust density versus the grain 
mass. 

of the electric field. 

The vertical dashed lines in Fig. |21 show the position of the sheath-presheath interface, 
z 8 , for 77 7^ 0, which is conditionally defined as the point where the ion flow velocity is equal 
to the ion sound velocity. The dependence of the coordinate z s on the dust parameters is 
rather weak. For this particular example, the dust layer levitates in the supersonic ion flow, 
zi < z s ] with reducing grain weight, /i, it moves up to the subsonic area, Z\ > z s . 

Fig-Elshows the dependence of the surface tension on the dust layer density, rj, for various 
values of the grain mass, \i. The curves 1-3 in Fig. |H] correspond to the dust layer levitating 
in the supersonic flow, z\ < z s , the curve 4 corresponds to lighter dust grains situated in 
the presheath area, zi > z s . The depicted curves show that if the dust layer is located in 
the sheath, the surface tension increases with r), while with the dust layer shifted to the 
presheath, the surface tension may become a decreasing function of rj. This is also shown in 
Fig. |H where the dependence of the derivative s / (^)| r?=0 on the grain mass, /i, is plotted. 

As we have already mentioned, the dependence of the surface tension on the dust density 
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determines the force acting upon the edge of the semi-infinite layer. If s(r]) > s(0) then the 
force f x = s(0) — s(rj) is negative, as shown in Fig. [H that is, the surface tension at the 
plasma-wall interface tends to expel the dust. Otherwise, the dust is drawn in the sheath. 
As it is seen from Figs. 13141 in dependence on the grain mass and the dust density both 
alternatives are possible. 

Here, we treat dust layer as a rigid incompressible medium. However, the short-range 
interaction between the dust grains result in the surface pressure, p s (cr), which also depends 
on the density [l^. The surface tension (fTTj) reduces the surface pressure, that is, the net 
equation of state of the dust layer is given by p = p s (o~) — s(a). Since p s (o~) exponentially 
tends to zero with decreasing a, the compressibility, dp/ da, of a rarefied dust layer is layer 
is negative if s'(0) > 0. Such a medium is unstable. The detailed investigation of this 
instability and arising structures will be discussed elsewhere. 

To summarize, charged dust may be treated as a surface-active substance that is capable 
of altering the surface tension at the plasma-electrode interface. Unlike usual surfactants, 
the dust increases the surface tension in a wide range ofparameters. Although we did not 
discuss here the emergence of the two-dimensional void [3j, it should be pointed out that a 
similar effect is well-known in liquids. To observe a void, it is enough to drop some liquid 
soap at the surface of water with floating fine powder. 

This study was supported in part by the Netherlands Organization for Scientific Research 
(grant no. NWO 047.008.013) and the Russian Foundation for Basic Research (project no. 
02-02-16439). 
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